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Abstract 

We consider a self-similar phase space with specific fractal dimension d being dis- 
tributed with spectrum function f(d). Related thermostatistics is shown to be 
governed by the Tsallis formalism of the non-extensive statistics, where the non- 
additivity parameter is equal to f{q) = l/r(q) > 1, and the multifractal function 
T (l) = Q.dq — f(d q ) is the specific heat determined with multifractal parameter 
q G [l,oo). In this way, the equipartition law is shown to take place. Optimiza- 
tion of the multifractal spectrum function f(d) derives the relation between the 
statistical weight and the system complexity. It is shown the statistical weight ex- 
ponent r(q) can be modeled by hyperbolic tangent deformed in accordance with 
both Tsallis and Kaniadakis exponentials to describe arbitrary multifractal phase 
space explicitly. The spectrum function f{d) is proved to increase monotonically 
from minimum value / = — latd = 0to maximum one / = 1 at d = 1 . At the same 
time, the number of monofractals increases with growth of the phase space volume 
at small dimensions d and falls down in the limit d — > 1 . 
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1 Introduction 



A generalization of the statistical mechanics onto the non-extensive thermo- 
statistics is known to be based on the deformation procedure of both logarithm 
and exponential functions [1,2,3]. The simplest way to introduce these func- 
tions into the thermostatistics scheme is to consider the equation of motion 
for dimensionless volume 7 = T/(2iih) 6N of the supported phase space (H, N 
being Dirac-Planck constant and particle number). In the course of evolution 
of the ensemble with statistical weight w — 10(7) and entropy S, the 

variation rate of the phase space volume is governed by the equation [3] 

d7 dS 

dt =w{l) W W 

Following from here relation dS = d r y/w{ r y) gives the entropy corresponding 
to the whole statistical weight W: 

y(W) 
7(1) 

Here, we take into account that entropy of a single state W = 1 vanishes, i.e., 
5(1) =0. 

In the case of the smooth phase space, one has trivial relation 111(7) = 7 whose 
insertion into Eq.(2) arrives at the Boltzmann entropy S = In IT. However, 
complex systems have fractal phase space with the dimension D < 6N, so that 
relation between the statistical weight and the corresponding volume should 
be generalized by the power-law dependence 

w(j) = 7 d (3) 

where the specific fractal dimension d = D/6N < 1 is introduced as the 
exponent. Insertion of Eq.(3) into the integral (2) gives the expression/ 1 "] 

S(W) = d\n 2 _j(W), \n q (x) = ^-^± (4) 

1 -Q 

which is reduced to the Tsallis logarithm ln q (x) where the non-additivity pa- 
rameter q is replaced by the difference 2 — d with d = 1/d > 1 being the 
inverse value of the specific fractal dimension d of the phase space. Naturally, 
this expression gives the Boltzmann entropy in the limit d — > 1. 

1 This expression is equivalent to Eq.(16) in Ref.[3] since W in our manuscript 
denotes M given by (11) in [3]. 
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Above formalism is based on the proposition that the phase space is related 
to a monofractal set determined by single dimension d. However, the consid- 
erations [4,5,6] show that a complex system behaviour can be determined by 
the phase space geometry, being much more complicated, in particular mul- 
tifractal. In this connection, we aim to generalize the Tsallis thermostatistics 
onto the multifractal phase space with a spectrum f(d). Such a generalization 
for arbitrary distribution f(d) is carried out in Section 2. Related discussion 
shows that physical representation of the thermostatistics based on the multi- 
fractal phase space demands of the passage from input distribution to escort 
one. An optimization procedure of the spectrum f(d) is considered in Sec- 
tion 3 to derive the relation between the statistical weight and the system 
complexity. In Section 4 we show that the monotonically increasing mass ex- 
ponent r(q), being free energy of the multifractal set [7], is presented by the 
hyperbolic tangent deformed in accordance with both Tsallis and Kaniadakis 
procedures, which allow for to describe explicitly arbitrary multifractal phase 
space. Section 5 is devoted to consideration of the multifractal spectrum f(d) 
which determines the number of monofractals within the multifractal with the 
specific dimension d. Section 6 concludes our consideration. 



2 Thermostatistics of multifractal phase space 



According to the self-similarity condition the specific statistical weight of the 
system under consideration is given by the power law function [8] 

wM = ^ (5) 

where q is the multifractal exponent, d = D/QN < 1 is the specific fractal 
dimension. This function should be multiplied by the number of monofractals 
with dimension d 

■Wr) = i~ m (6) 

which are contained in the multifractal whose spectrum is determined by a 
function f(d). As a result, whole statistical weight, being the multifractal 
measure, takes the form 

i i 
w,( 7 ) = J w q ( 1 )Af d ( 1 )p(d)dd = J ^ d -^p(d)dd (7) 



where p(d) is a density distribution over dimensions d. Using the method of 
the steepest descent, we arrive at the power law 

w q (>y) ~ Y (q) (8) 
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which generalizes the simplest relation (3) due to replacement of the bare 
fractal dimension d by the multifractal function 



r(q)=qd q -f(d q ), (9) 

being the mass exponent [8]. Here, the specific fractal dimension d q relates to 
given parameter q to be defined by the following conditions of the steepest 
descent method: 



dd 



-q ^ 

d=d a 



<o. (10) 

d=d a 



Above consideration shows that the passage from the monofractal phase space 
to the multifractal one is obtained by replacement of the single dimension d by 
the monotonically increasing function r(q), such as r(0) = —1 and r(l) = 0. 
Limit behaviour of the function r(q) is characterized by the asymptotics [8] 

ra(g-l) at 0<g-l<l, r~l at g^oo. (11) 

A physical domain of the q parameter variation is bounded by the condition 
q > 1 which ensures positive values of the function < r(g) < 1 to guarantee 
growth of the specific statistical weight (8) with increasing the phase space 
volume. 

According to the entropy expression (4) we can use well-known Tsallis for- 
malism of the non-extensive statistical physics where the difference 2 — f (q) 
with f(q) = l/r(q) > 1 plays a role of the non-additivity parameter. Thus, 
the entropy in dependence of the probability distribution Pj has the form [1] 

Wg pf-l _ , 

5 * = -E p H^ri- ^ 
i=i t i 

where the statistical weight W q is related to given value q of the multifractal 
exponent. With accounting the normalization conditions and the definition of 
the internal energy E q 

Wg Wg 

£p = i, £ ff = 5>i?<«\ (13) 
1=1 1=1 

the expression (12) arrives at the generalized distribution over energy levels 
Si as follows: 



Pi = Z q 1 expf(q) (~P £ i) . Z <1 = E eX Pr( 9 ) (~P £ i) ■ ( 14 ) 

i=l 

Here, Z q is the partition function, f3 is Lagrange multiplier, not being the 
physical temperature, and the deformed exponential function is determined 
by the expression 
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exp f {x) = < 



[1 + (f - l)x] r-i at 1 + (f - l)x > 0, 

(15) 

otherwise. 



Parameter q characterizes here the multifractal spectrum through the function 
r(q) and should no be confused with the non-additivity parameter of Tsallis 
thermostatistics, which is denoted here by f(q) = l/r(q). 

Thermodynamic functions of the model under consideration can be found 
according to the Tsallis non-extensive scheme [1] . However, related expressions 
are very cumbersome even in the simplest case of the ideal gas [9,10,11] and 
take the usual form only within the slightly non-extensive limit [12]. At the 
same time, developed scheme allows for to use thermodynamic formalism of 
multifractal objects [13], within the which the multifractal exponent q plays a 
role of a state parameter. If the dependence r(q) has some singularities, then 
variation in q may arrive at phase transitions. It is worthwhile to stress that 
developed scheme arrives directly at related singularities of thermodynamic 
functions type of the internal energy (see below Eq.(23)), the entropy (cf. 
Eq.(4)) 

S q = f(q) ln 2 _ f((?) (W q ) , f(q) = l/r(q) (16) 
and the free energy 

F q = E q - TS q . (17) 

According to Ref. [9], the physical distribution is not the input probability 
(14), but the escort one 

V(si) = -4r ■ (18) 

It corresponds to the condition 

w q 

J2(e t -E q )P f ^(e l ) = (19) 

i=l 

instead of the second equation (13). In difference of the distribution (14) re- 
lated probability 

Vie,) = Z' 1 ex Pr(9) l-f(q) (e t - E q ) /T] (20) 
is determined with the physical temperature T. 

In the case of continuous energy spectrum characterized with the density dis- 
tribution p(e), the internal energy related to the condition (19) takes the form 

oo 

E q = J eV(e)p(e)de. (21) 
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Extreme value of E q is reached at the condition 



^M~-^M (22) 

where prime denotes differentiation over e. Usually, the density function is 
reduced to the power law p(e) ~ e cN , c ~ f , so that p (e)/p(e) ~ cN/e. Then, 
with using the distribution (20), the condition (22) taken at e = E q arrives at 
the equipartition law 

E q = cr(q)NT (23) 
where the value cr(q) is the specific heat. 



3 Optimization of multifractal spectrum 



Up to now, we suppose that the multifractal spectrum f(d) is arbitrary. If it 
is optimized at normalization condition 



/ f(d)dd 



I- 



(24) 



one has to minimize the expression 

7(W,) 



S q {f(d)} = 



7(1) 



7 9d - /(d) p(d)dd 



-l 



d7 



f(d)dd-l 



(25) 



where we take into account Eqs. (2), (7), and E determines the Lagrange 
multiplier. As a result, we arrive at the equality 



i(w q ) 



J 7 {[?d-/(d)]-2r( ? )} ln7 ^ 



5! 



(26) 



7(1) 



whose integration gives, with accounting Eq.(8), the transcendental equation 



\ lZr(q)T q (d)] 2 - [T q (d) ln(W q ) - 1] - 1 = 0, (27) 

T q (d) = f(q) {1 - 2r(q) + [qd - f(d)]} . (28) 

This equation is written in the form, when can be used either the spectrum 
function f(d) or the exponent dependence r(g). In the latter case, we find 
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initially the dependence q(d) from the equation 

= d, 



dr 
dq 



(29) 



q=q(d) 



being conjugated to Eq.(lO). Then, inserting this dependence into Eq.(28), we 
arrive at the trivial expression 



T(q) = T q {d{q)) = f{q) 



(30) 



whose using in Eq. (27) allows for to determine the dependence of the statistical 
weight W q versus the complexity S at given function r(q). A typical form of 
this dependence at the mass exponent 



q-1 



(31) 



is shown in Fig.l. It is seen the statistical weight increases monotonically as 
a b 




Fig. 1. Dependencies of the statistical weight on the complexity (a) and multifractal 
exponent (b) at the exponent (31) (numbers near curves on the left panel show values 
q, on the right one - values S) 

with growth of the complexity, so with increasing multifractal exponent. 

In the limit of smooth phase space, when q — > oo, d — > 1, T(q) — > 0, one 
obtains the usual expression for the statistical weight of the complex system 



(32) 



which is determined by the specific complexity per one particle. At small 
deviation off the minimum complexity (E — <^ T,^) and light multifrac- 
tality (1 — r(q) <^ S^ 1 ), linearized equation (27) gives 



W q ~ 



W T{q) { 1 



r(q) 

W T (q) = exp |r(g)Soo 



l--(l-r( g ))E 00 
1- hl-r( g ))E 00 



(S-Soo)}, 

| , l«g<oo. 



(33) 
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In the opposite case r(q) <C 1, one has within logarithmic accuracy 



ar(q)J: 2 , a ~ 1, 9 - 1 < 1. (34) 



As show above findings, optimization of the multifractal spectrum, obeying 
the normalization condition (24), gives the dependence of the statistical weight 
W g versus the system complexity £ at given multifractal exponent r(q). Nat- 
urally, relations (27) and (30) allows for to solve the inverse problem - to 
find the dependence r(q) at given function W q (Yi). However, definition of the 
dependence W g (£) leads to very complicated problem. It is more convenient 
to use a modeling function r(q) bounded with asymptotics (11) and then to 
find the statistical weight W q . Within this algorithm, in the following Section 
we model the multifractal spectrum on the basis of the procedure of both 
Tsallis and Kaniadakis deformations. It is appeared, such deformations give 
the whole set of functions f(d) to present all possible types of multifractal 
spectra. 



4 Analytical modeling multifractal spectrum 



In the simplest case, we can take the function r{q) in the form 

r = tanh [V(q - 1)] (35) 

being determined by parameter V > and argument q G [1, oo]. According to 
Fig. 2 related multifractal dimension function [8] 

D q = (36) 

monotonically decreases from the maximum value Dq = T> to the minimum 
one -Dqo = with increase in q. However, the maximum value of the fractal 
dimension D q is fixed by the magnitude D = 1, so that one should put V = 1 
in the dependence (35). As a result, it takes quite trivial form. 

Due to the r(g) function increases monotonically within narrow interval [0, 1], 
one has a scanty choice of its analytical models. To set a possible represen- 
tation of r(q) one can use a deformation of the hyperbolic tangent (35) at 
T> — 1. By now, two analytical procedures of such deformation are extensively 
popularized. The first of them is based on the Tsallis exponential form [1] 



exp K (x) = < 



(1 + kx) 1,k at 1 + kx > 0, 

V (37) 
otherwise 
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100 



Fig. 2. Multifractal dimension function related to the dependence (35) (curves 1, 2, 
3 correspond to V = 0.1, 1.0, 10). 

where deformation parameter n takes positive values. The second procedure 
has been proposed by Kaniadakis [14] to determine the deformed exponential 
form 



With using these definitions, the deformed tangent (35) takes the form 



where the multifractal exponent q varies within the domain [1, oo]. 

The g-dependencies of the multifractal exponent r(g) and its inverse value 
f(q) = 1/V(g) are shown in Fig. 3 at different magnitudes of both Tsallis and 
Kaniadakis deformation parameters k. (The k values are picked out in such 
a manner to cover uniformly the panels of Figs. 3-5 with related curves.) It 
is principally important, the first of these deformations arrives at more fast 
variations of both exponents r(g) and f(q) in comparison with non-deformed 
hyperbolic tangent tq = tanho(g — 1), whereas the Kaniadakis deformation 
slows down these variations with n increase. 

A characteristic peculiarity of the Tsallis deformation consists in breaking 
dependencies r(q), f(q) in the point go — (1 + k)/k where the second terms 
in both numerator and denominator of the definition (39) take the zero value. 
As a result, the multifractal exponent (9) takes the asymptotics 




(38) 



T K (q) = tanh K (g - 1) 



exp„(g - 1) - exp K (l - q) 
exp K (g- 1) +exp K (l - q) 



(39) 
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3.0 




Fig. 3. The (/-dependencies of the multifractal exponent r and its inverse value 
f = 1/r (solid line corresponds to k = 0, dashed curves relate to the Tsallis 
deformation with k = 0.7,1; dotted lines correspond to the Kaniadakis one at 
k = 0.7, 1,3,10). 



r, 



(Ts) 



' (g-l)-^(g-l) 3 
' l- 2(«/2) 1/K (g -g) 1/K 



at < g - K 1, 
at < q - q < go- 



(40) 



For k = 1 the dependence r} SJ (g) takes the simplest form: Ti' — q — 1 
at 1 < g < 2, and r^ 5 - 1 = 1 at g > 2. It is worthwhile to note, the Tsallis 
deformation parameter can not take values k > 1 because these relate to the 
fractal dimensions D q > 1 at g 7^ 0. 

In the case of the Kaniadakis deformation, the multifractal exponent r K (g) 
varies smoothly to be characterizes by the following asymptotics: 



(g-l)-^(g-l) 3 at < g - 1 « 1, 
1 - 2 [2/c(g - l)p 2/K at «(g - 1) > 1. 

In contrast to the Tsallis case, here the deformation parameter can take arbi- 
trary values to give the simplest dependence (31) at k — 1. 

The fractal dimension (36) as a function of the q exponent falls down mono- 
tonically as shown in Fig.4a. According to Eqs.(40), in the case of the Tsallis 
deformation, one has a broken dependence D(q), being characterized by the 
asymptotics 



T {K) ^ J 

K 
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Fig. 4. Spectra of fractal (a) and specific (b) dimensions of the phase space (solid 
line corresponds to k = 0, dashed curves relate to the Tsallis deformation with 
k = 0.7, 1; dotted lines correspond to the Kaniadakis one at k = 0.7, 1,3, 10). 



D (Ts) 



1 - ^ (q - 1 
li"2(«/2) (/ , 



i/k (qo-<?) 1A 



at < q - K 1, 
at < g - q < <?o- 



(42) 



In the limit case k = 1, the phase space is smooth (D q Ts ^ = 1) within the 
interval 1 < q < 2. For the Kaniadakis deformation, the fractal dimension 
is given by smoothly falling down curve whose slope increases with the 
deformation parameter growth. According to Eqs.(41), in this case, one has 
the asymptotics 



q ~ 



l-^(q-l) 



at < q - K 1, 
(q - l)' 1 - 2 (2kY 2/k (q - l)-(-+ 2 )/« at «(g - 1) > 1. 



(43) 



At k — 1, the typical dependence related to Eq.(31) is of the form 



D (K) 



1 + (?-!) 



-1/2 



(44) 



5 Multifractal spectrum 



At given multifractal exponent r(g), the spectrum function f(d) is denned by 
the Legendre transformation (9) where the specific dimension reads 

d , = |. (45) 

As shows Fig. 4b, the dependence d q has monotonically falling down form to 
take the value d q = at q > q = (1 + /c)//c for the Tsallis deformation. In 
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this case, asymptotical behaviour is characterized by Eqs.(40), according to 
which one obtains 



1 - (1 - k 2 ) (q - i f at 0<g-Kl, 

(K/2f~ K)/K (g -qf-«V« at < q - q « q . 



(46) 



In the limit n — > 1, the dependence d <yTs \q) takes the step-like form being 
d q — 1 within the interval 1 < q < 2 and d ? = otherwise. 

For the Kaniadakis deformation, Eqs.(41) arrive at the asymptotics 



1 - 



2 + K Z 



(9-1) 



-3/2 



2 2(«-l)/« !)]-(2+«)/« 

The typical behaviour is presented by the dependence 

df ] = [l + (q - I) 2 
related to k — 1. 

The multifractal spectrum is defined by the equality 

/(d) = dq d - r(q d ) 



at < q - 1< 1, 
at «(g-l)»l. 



(47) 



(48) 



(49) 



being conjugated to Eq.(9). Here, the specific multifractal exponent q d is de- 
termined by Eq.(45) which arrives at the limit relations (46), (47). With their 
using, one obtains the asymptotics 



./ 



(Ts) 



d-\(l-K 2 )- l l 2 {l-dfl 2 



at < 1 - d< 1, 



- 1 + 2(1/k- l)d l ^ l -^\ + (1 + l/n)d at d<l 
for the Tsallis deformation, and the relations 



(50) 



./ 



(K) 



d 



1 



-1/2 



1 - cf) 3 / 2 at < 1 - d < 1, 



3 V 2 , v- - - ^ (51) 

1 _ 2 (*-4)/(2+«)( 1 + 2//t)d 2 /( 2+K ) + d at d < 1, 



characterizing the Kaniadakis deformation. 

As shows Fig.5, for finite deformation parameters k < oo, a spectrum function 
increases monotonically, taking the minimum value / = — 1 at d = and the 
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1.0 




/ 



0.5 



0.0 



-0.5 



-1.0 



Fig. 5. Spectrum function of the multifractal phase space (solid line corresponds to 
k = 0, dashed curves relate to the Tsallis deformation with k = 0.7, 1; dotted lines 
correspond to the Kaniadakis one atre=l,3,10,oo). 

maximum one / = 1 at d — 1. Besides, the derivative /' = df/dd equals to 
f'{0) = oo on the left boundary and /'(I) = 1 on the right one. It is significant, 
the whole set of the spectrum functions is bounded by the limit dependencies 
j?(Ts) — 2d — 1 and = d, the first of which relates to limit magnitude of 
the Tsallis deformation parameter k — 1 whereas the second one corresponds 
to the Kaniadakis limit k — oo. A typical form of the spectrum function is 
presented by the dependencies 



It may seem, at the first glance, that negative values of the spectrum function 
f(d) has not a physical meaning. To clear up this problem, let us take the set of 
monofractals with the specific dimension d = 0. Obviously, such monofractals 
relate to the whole set of the phase space points, whose number equals to the 
dimensionless volume 7. Just such result gives the definition (2) in the point 
d — where / = — 1. On the other hand, in opposite case d — 1 where / = 1, 
we obtain the very thing the number of monofractals with volume 7 equals 
A/1 = 7 _1 to give one multifractal of the same volume 7. At the same time, a 
single monofractal is contained in the multifractal at condition f(d) = which 
takes place at the specific dimension do whose dependence on the deformation 
parameter k is shown in Fig. 6. 

The dependence of the number M of monofractals containing in the phase 
space volume 7 related to the multifractal with the specific dimension d is 
shown in Fig. 7. It is seen, the number M increases with the 7 volume growth 



f(K) = 




(52) 
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0.2 0.4 0.6 0.8 1.0 




10- 2 10 1 10° 10 1 10 2 10 3 10 4 

K 

Fig. 6. Dimension do related to the condition f(d) = as function of the parameter 
k (curve 1 corresponds to the Tsallis deformation, curve 2 - to the Kaniadakis one; 
positive values f(d) relate to the domain d > do). 




Fig. 7. The number Af of monofractals within the phase space volume 7 related 
to the specific dimension d at the deformation parameter k = 10 (different levels 
Af = const are shown with thin curves, thick lines relate to Af = 1). 



at small dimensions d, whereas in the limit d — > 1 the dependence A/*(7) 
becomes falling down to give infinitely increasing numbers Af at 7 — * 0. The 
speed of such increase growths monotonically with both decrease of the Tsallis 
deformation parameter k and increase of the Kaniadakis one. 
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6 Conclusions 



As above consideration shows, the statistical mechanics of self-similar complex 
systems with phase space, whose specific fractal dimension d is distributed 
with spectrum f(d), is governed by the Tsallis formalism of the non-extensive 
thermostatistics. In this way, the role of non-additivity parameter plays inverse 
value of the multifractal function r(q) = qd q — f(d q ) which monotonically 
increases, taking value r = at q = 1 and r ~ 1 at q — > oo (the latter limit 
relates to the smooth phase space, where r = 1). The multifractal function r(q) 
is reduced to the specific heat to determine, together with the inverse value 
f(q) = l/r(q) > 1, both statistical distributions and thermodynamic functions 
of the system under consideration. At given function r(q), optimization of 
the normalized multifractal spectrum f(d) arrives at the dependence of the 
statistical weight on the system complexity. 

It is shown the whole set of monofractals within a multifractal related to the 
phase space, which gives the support of a generalized thermostatistics, is mod- 
eled by the mass exponent r(q) that determines the statistical weight (8) at 
given volume 7. To be the entropy (2) concave, Lesche stable et cetera, the ex- 
ponent r(q) should be a function, monotonically increasing within the interval 
[0, 1] at multifractal exponent variation within the domain [1, 00]. The simplest 
case of such a function gives the hyperbolic tangent r = tanh(g — 1) whose 
deformation (39) defined in accordance with both Tsallis and Kaniadakis ex- 
ponentials (37), (38) allows for to describe arbitrary multifractal phase space 
explicitly. At the same time, the Tsallis deformation arrives at more fast varia- 
tions of the statistical weight exponent r(q) in comparison with non-deformed 
hyperbolic tangent, whereas the Kaniadakis one slows down these variations 
with increasing the deformation parameter k. All possible dependencies r(q) 
are bounded from above by the linear function r^ Ts ^ = q— 1 at q E [1,2] which 
is transformed into the constant r = 1 at q > 2. This dependence relates to 
the smooth phase space within the Tsallis interval q G [1, 2]. 

The dependence (6) of the number of monofractals within the phase space 
volume 7 related to the multifractal with the specific dimension d is determined 
by the spectrum function f(d). This function increases monotonically, taking 
the minimum value f — —1 at d — and the maximum one / = 1 at d — 1; 
besides, its derivative equals /'(0) = 00 on the left boundary and /'(l) = 1 
on the right one. The whole set of the spectrum functions is bounded by the 
limit dependencies f^ Ts ^ = 2d — 1 and f( K ) = d, the first of which relates to 
limit magnitude of the Tsallis deformation parameter k — 1 and the second 
one corresponds to the Kaniadakis limit k — 00. The number of monofractals 
within the multifractal increases with the 7 volume growth at small dimensions 
d and falls down in the limits d — > 1 to give infinitely increasing at 7 — > 0. 
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